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Abstract: The methods of macroscopic averaging used to derive the macroscopic Maxwell 
equations from electron theory are methodologically incorrect and lead in some cases to a 
substantial error. For instance, these methods do not take into account the existence of a 
macroscopic electromagnetic field E#, H# generated by carriers of electric charge moving in 
a thin layer adjacent to the boundary of the physical region containing these carriers. If 
this boundary is impenetrable for charged particles, then in its immediate vicinity all carriers 
are accelerated towards the inside of the region. The existence of the privileged direction of 
acceleration results in the generation of the macroscopic field E#, H#. The contributions to 
this field from individual accelerated particles are described with a sufficient accuracy by the 
Lienard-Wiechert formulas. In some cases the intensity of the field E#, H# is significant not 
only for deuteron plasma prepared for a controlled thermonuclear fusion reaction but also for 
electron plasma in conductors at room temperatures. The corrected procedures of macroscopic 
averaging will induce some changes in the present form of plasma dynamics equations. The 
modified equations will help to design improved systems of plasma confinement. 
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1. INTRODUCTION 

The methods of macroscopic averaging used to derive the macroscopic Maxwell equations 
are based on two false mathematical assumptions and lead sometimes to false results. This will 
be explained in the present section. 

It has been directly shown in my earlier paper (Szalek 1997) that the correct procedure of 
macroscopic averaging of the electromagnetic field generated by particles from a small macro- 
scopic region P produces a field not always satisfying the macroscopic Maxwell equations. This 
follows from a simple example discussed in that paper or from comment 2 on p. 100 there. 
This result may be considered a consequence of the fact that each contribution to the field 
from an individual carrier leaving or entering P is discontinuous in time. It is different from 
zero when the generating carrier is in P (in the retarded time) and is equal to zero otherwise. 



Therefore, none of these contributions satisfies the pertinent microscopic Maxwell equations, 
which invalidates, e.g., the formulas (45), (46) on p. 67 and (84), (85) on p. 113 in the book of 
de Groot (1969) as well as formulas (8) on p. 24 and (62) on p. 256 in the book of de Groot 
and Suttorp (1972). These formulas average the microscopic Maxwell equations with pseudo- 
solutions introduced to these equations, where the pseudo-solutions are the contributions to 
the electromagnetic field from individual carriers. The aim was to obtain the macroscopic elec- 
tromagnetic field being the solution of the macroscopic Maxwell equations. This procedure is 
erroneous as the pseudo-solutions do not satisfy the pertinent microscopic Maxwell equations 
if the boundary of P is penetrable for carriers. 

If a proof is erroneous, it need not necessarily mean that its thesis is incorrect. That, 
however, is not the case here, and the proper solutions of the macroscopic Maxwell equations 
not always describe correctly the macroscopic electromagnetic field. The differences between 
these solutions and the real macroscopic fields are sometimes surprisingly great. The proof is 
given in this paper in some cases when the macroscopic electric current does not flow through the 
boundary of the region P. Such an outcome has been already suggested by results of my earlier 
paper concerning some cases when the macroscopic current flows through the boundary of P. 
The results obtained in the present paper make it possible to propose in Section 2 a convenient 
procedure which will allow to calculate the correct value of the averaged electromagnetic field 
acting on charged particles in material media. 

The derivation of the macroscopic Maxwell equations is also accompanied by another 
methodological error. Namely, the pertinent considerations take into account the distribu- 
tion functions of particles depending only on space coordinates and on the velocities of carriers 
of electric charge. They disregard the distribution of carrier accelerations. Consequently, the 
particles with the same velocities but different accelerations are treated identically, although 
the fields generated by these particles are different. Such an approach may be justified when 
in a small macroscopic region, or in a thin layer adjacent to a privileged surface, the average 
acceleration of all carriers is equal to zero. This approach is incorrect when we consider a 
thin layer adjacent to the impenetrable boundary of a physical medium, where the average 
acceleration is different from zero and is directed towards the inside of the medium. It is also 
erroneous when we take into account an macroscopically inhomogeneous medium. 

The considerations of this paper are based on the existence of the macroscopic electro- 
magnetic field generated by carriers of electric charge accelerated in the immediate vicinity 
of impenetrable boundaries. This field has not been taken into account in the pertinent ar- 
gumentation of other authors. Its intensity is sometimes astonishingly great, which indicates 
that the whole problem of the macroscopic electromagnetic field in material media should be 
meticulously reconsidered. 

This paper is an abbreviated and simplified version of my earlier paper "Plasma dynamics 
and boundary phenomena in macroscopic averaging" formulated in 2002/2003. In that paper 
results concerning the correct calculation of the macroscopic electromagnetic field have been 
obtained in an exact and systematic manner in the general case of curved boundaries and real 
physical media. Certain considerations there, however, seem rather difficult for some physicists, 
and the exactness of the argumentation has made that paper relatively long (25 pages). To 
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make the argumentation of the present paper shorter and simpler, we idealize to some extent 
the properties of the media in the vicinity of the impenetrable boundary. We also assume that 
the boundaries are plane. Nevertheless, the results obtained here agree with the results of that 
earlier paper. These of the more important results of the latter which have not been derived 
here, are repeated without proofs in Section 5. We have also retained some methodological re- 
marks which are given in Section 6. That earlier paper is available from the author on request 
electronically in a form of a pdf file. 

2. INITIAL CONSIDERATIONS 

Consider a macroscopic bounded region P containing a plasma or a typical metal without 
significant magnetic properties. We denote the boundary of P by dP and assume that dP is 
composed of two sufficiently regular surfaces d\P and d 2 P ■ The surface d 2 P is impenetrable 
for carriers of electric charge. 

Let J and p be the macroscopic densities of, respectively, electric current and charge in 
P, and let the normal to dP component of J vanish identically on dP. Let E , H denote 
respectively the electric and magnetic field being the solution of the Maxwell equations with 
the densities J and p, and satisfying the proper boundary conditions. If J = 0, p = in P, then 
E = 0, H = 0. According to the present opinions the field E , H is equal to the macroscopic 
electromagnetic field E, H generated by particles from the region P: 

E = E , H = H . (1) 

The field E, H is the sum of contributions to the macroscopic electromagnetic field from all 
particles in P. We shall show that in order to describe correctly this field, one should take 
into account in some cases not only the field E , H , but also two other macroscopic fields 
which will be denoted by E B , H B and E T , H T . The field E B , H B is generated by carriers 
accelerated in the immediate vicinity of the boundary d 2 P ■ For all carriers this acceleration is 
directed towards the interior of P and is induced by forces which make d 2 P impenetrable. The 
existence of the privileged direction of acceleration causes that the sum of contributions to the 
field E#, Hb from particular carriers does not vanish. It will be shown that the intensity of the 
field E B , H B can be significant even in some cases of electron plasma in conductors at room 
temperatures. We shall refer to this field as "the boundary acceleration electromagnetic field." 

If E = 0, H = 0, then the field E^, is the sum of contributions, respectively, to the 
macroscopic electric and magnetic field, generated by all particles from P with the exception of 
the contributions used to calculate E^, H#. If Eo 7^ or Ho 7^ 0, then we define that this sum 
is equal to E + E^, H + Hy We shall call the field Et, Ht "the temperature electromagnetic 
field." 

Let [E, H] be an ordered pair and [E, H] = [e, h] if and only if E = e and H = h. Let 

G = [E,H], Go = [E ,H ], Gb = [Eb,Hb], Gt = [E t ,Ht}- (2) 
Using this notation we may write the correct expression for the field G in the form 

G = Go + Gb + Gt- (3) 
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This and the next equations concern the fields G, Gb, Gt in observation points lying outside 
the region P and its immediate vicinity. The fields Gb and G^ may be ignored when 

G B + G T = 0. (4) 

We shall show, however, that in some cases equation (J1J) is not satisfied. 

First, consider the case when the plasma in P is homogeneous and isotropic. Let |G| = 
|E| + |H| where |E|, |H| are the lengths of E and H. It follows from the results obtained in the 
present paper that if d^P = dP, then \Gb\ vanishes or is relatively small. We may expect the 
same from |Gy| taking into account equation (J3J, the validity of which is in this case suggested, 
e.g., by the experimental fact that the electromagnetic field produced by carriers from any 
region P with impenetrable boundaries, and averaged in a short interval of time, is at the 
outside of P and its immediate vicinity equal to zero with a good accuracy for media under 
consideration if the macroscopic densities J and p vanish in P. 

The situation changes substantially when d?P 7^ dP and d?P is not a closed surface. In such 
cases \Gt\ remains roughly the same while (G^l can become surprisingly large. For instance, 
one impenetrable face of a copper cube with the length of edge equal to 1 cm produces at a 
temperature of 300 K the fields E#, which in some points at a distance of 1 m from the 
center of the cube have the intensity of about, respectively, 10 11 Vm" 1 and 10 8 Am _1 . If the 
cube were filled with a rarefied deuteron plasma prepared for a controlled thermonuclear fusion 
reaction, then the corresponding intensities would be about 10 6 Vm -1 and 10 4 Am _1 . We can 
obtain then from (jSJ) 

|G-G | = |G B + G T | >0. (5) 

In particular, we can get |G| 3> when J = 0, p = in P. A proof of inequality (jSJ) for Go = 
is given in Section 3. 

Inequality (jSJ) is interesting from the theoretical point of view, because it evidently disagrees 
with equations (0) and with some universally accepted opinions concerning the macroscopic 
electromagnetic field generated by large sets of charged particles in chaotic motion. There 
are two main reasons of the disagreement between the results obtained in the present paper 
and the pertinent results of other authors. First, the methods of macroscopic averaging of 
electromagnetic fields applied in the scientific literature do not take into account the existence 
of the macroscopic boundary acceleration field Gb generated in a locally microscopic layer 
adjacent to c^P. Second, it has not been noticed that for inhomogeneous or anisotropic media 
the temperature field Gt may be different from zero and attain considerable intensities. An 
additional reason of the disagreement is the discontinuity in time of the contributions to the 
macroscopic electromagnetic field from individual carriers, which has been already explained 
in Section 1. 

If the medium in P is inhomogeneous and/or anisotropic, then \Gb\ can become relatively 
very large even if d^P = dP or 82 P is a closed surface. This may have practical consequences. 
It will be explained in the following example. 

Suppose that we want to calculate an averaged electromagnetic field G' = [E', H'] acting on 
electrons and ions at a point O of a region P' with a boundary dP'. There are various methods 
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of calculating this field. All of them need some corrections because they do not allow for the 
existence of the fields Gb, Gt- We shall point out a method where the manner of correction 
is obvious. 

Let the medium in P' be inhomogeneous and the surface dP' be impenetrable for carriers. 
To facilitate the considerations we assume that the medium inhomogeneities are stable, as in 
inhomogeneous alloys of metals. Consider in P' a small subregion P$ with a boundary dPs 
such that O G Ps and lies possibly far from dPs- We may write down 

G' = G s + G (6) 

where the field Gs is generated by particles from inside the region P$, and the field G is 
produced by particles from inside the region P defined by P = P' — Ps- The field Gs may be 
calculated by means of summing the contributions from all particles in Ps with the exception 
of a privileged particle, while the field G is calculated by means of the macroscopic densities J 
and p in P. The standard practice is to assume that the latter field satisfies (JTJ), i.e. G = Go- 
On the ground of results obtained in this paper it can be shown that the assumption (JTJ) is 
not always correct even if Go = 0. This is a consequence of the facts that dP = dP' + dPs, 
and that the surface dPs is penetrable for carriers and does not generate a contribution to the 
field Gb. Let the contribution generated at dP' to the field Gb be G' B and let the normal 
to dPs component of J vanish on dPs- We denote by Gbs the contribution to the boundary 
acceleration field which would be produced at dPs by particles from P if dPs were impenetrable 
for carriers. In such a case the whole boundary dP of P would be impenetrable for carriers. 
Therefore, we obtain from (0J) 

G T + G' B + G BS = 0- (7) 
We receive from (J3J) when dPs is penetrable: 

G = G + G' B + G T . (8) 

From (J7j) and (jSJ) we get 

G = Go - Gbs. (9) 

An example is given in Section 4, which indicates that in inhomogeneous plasma one can obtain 
| Gbs | ^ as in (jSOj) . Consequently, instead of the incorrect formula (JTJ), we should use equation 
(JSJ). From (JHJ) and Q we get the correct expression for the field G': 

G' = G s + G - Gbs- (10) 

At this stage of the research one cannot exclude the possibility that the presence of the term 
Gbs i n can have a noticeable influence on the value of G'. Such an influence should cause 
some changes in the present form of plasma dynamics equations. One can expect the greatest 
values of |Gbs| in the regions where the strongest inhomogeneities are present, e.g., on a shock 
wave. 

In general, the field Gbs w ih depend on the derivatives of plasma density and temperature. 
Hence, its existence could cause some changes of values of the coefficients of partial derivatives 
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in the plasma dynamics equations. In particular, some coefficients equal to zero in the present 
form of these equations might become different from zero. Even small differences may affect, 
e.g., results concerning the growth and propagation of plasma instabilities. In order to examine 
the consequences of equation (fTU|) one should derive the explicit dependence of the field G bs on 
the macroscopic parameters of inhomogeneous plasma, and investigate relations between Gbs 
and the coefficients in the plasma dynamics equations. 

The proper methodology pointed out in this paper will be useful not only when one examines 
the consequences of equation (fTUj) but also when it is possible to investigate the influence of some 
additional forces generated by electrically charged particles. These forces have been discovered 
by the present author and are briefly described in an earlier paper (Szalek 1997). Allowing 
for all forces produced by electrically charged particles, one will be able to obtain, e.g., an 
adequate description of phenomena occurring in boundary regions of globular lightnings. Such 
a description will be connected with further changes in the plasma dynamics equations. The 
modified equations will help to design improved systems of plasma confinement. 

The temperature and boundary acceleration fields are generated in various media and should 
be taken into account in the completed macroscopic theory of electromagnetic fields based on 
the microscopic properties of electrons, ions, and atoms. This theory came into existence a 
century ago, although some of its important premises were known much earlier - compare, e.g., 
the early theories of dielectrics or Ampere's hypothesis that magnetism is due to micro-currents. 
First results concerning the derivation of macroscopic Maxwell equations from electron theory 
were obtained by Hendrik A. Lorentz (see, e.g., Lorentz 1902, 1904, 1915). A list of authors 
working on this problem and a historical survey were given by de Groot (1969). First results 
concerning the fields Gb and Gt were published by the present author (Szalek 1997). That 
paper contains, e.g., an estimate of the ratio |H B + H r |/|H | (using the nomenclature intro- 
duced here) in some cases when P is a cylinder and d^P is the lateral surface of this cylinder. 
Because of lack of funds for the pertinent research, my next paper on this subject was written 
with a delay of more than 5 years. 



3.1. The field of a moving charge 

We shall calculate in sections 3.2 and 3.3 the contributions to the macroscopic electromag- 
netic field from particular carriers using some approximate expressions resulting from the exact 
Lienard-Wiechert formulas (Lienard 1898, Wiechert 1900). In the SI systems of units these 
expressions may be written down in the form 



where e(x, t), h(x, t) are, respectively, the electric and magnetic field generated by a moving 
point charge Q and observed at a point x at time t, 



3. A PROOF OF INEQUALITY © 









(12) 
(13) 
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where Ri(t) = x — r(t); x, r are, respectively, the position vectors of the observation point 
and the point charge in a Cartesian rectangular coordinate system X with the coordinates xi, 
X2, X3', -Ri = |Ri|, ni(i) = Ri/i?i; a, u are, respectively, the acceleration and velocity of 
the point charge; e = 8.85 x 10~ 12 F m _1 , c is the speed of light, and 

U = g(47T£ c 2 )- 1 , w = {l-n 1 - u/c)- 1 . (14) 

We have neglected in (|T2"j) and (fT3*Jl the time retardation, and in (JT2*Jl we have retained only 
the first term of the expansion with respect to the powers of u/c. The exact formulas may be 
found in various books and monographs (see, e.g., Clemmow and Dougherty 1969). They are 
also used in my previous pertinent paper. 

3.2. Introductory considerations 

Consider two identical macroscopic cubes P^ (i = 1, 2) containing a suitable homogeneous 
medium with the macroscopic properties not depending on time. We assume for the sake of 
simplicity that the ions of the medium are at rest in pW ( se e, however, remark 2 in Section 5). 
The cube P^) is inside the medium and its boundary dP^ = diP^ is fully penetrable for 
electrons. The cube P^ is adjacent to the boundary of the medium and one of its faces 
(denoted by c^P'- 2 -') is impenetrable for electrons while the remaining boundary diP^ of P^ 
is penetrable. 

Let Gg , Gj) denote, respectively, the macroscopic electromagnetic fields, the bound- 
ary acceleration electromagnetic fields, and the temperature electromagnetic fields generated 
by charged particles from pM, and let Gq be the proper solutions of the macroscopic Maxwell 
equations in pW (i = 1,2). These fields will be calculated in observation points lying outside 
the cubes PW, sufficiently far from their boundaries. We shall prove that at least one of the 
fields satisfies the formula (J5J) 

Assume that the macroscopic current and charge densities in P^ are identically equal to 
zero. Therefore, 

G^=0, z = l,2. (15) 

Let e a , h a , e u , h u be given by (fT2j) and (JT3J). The equation (fT5j) being satisfied, we can calculate 
the field G? by summing up the contributions [e a ,h a ], [e u ,h u ] from all particles in P^ (i.e. 
from electron carriers and ions in PW, the contributions from ions being given in our case by 
the Coulomb field). These contributions are averaged in a short interval of time when necessary. 

Assume now that the barrier of forces making the face c^P^ 2 -' impenetrable is sufficiently 
strong and steep to decelerate and accelerate the carriers within a distance k which is small 
in comparison with the average distance between carriers in P^ . Let AP( 2 ) be a boundary 
layer with a thickness k adjacent to the face c^P 1 - 2 ), and AP^) be a corresponding layer in 
P^ l \ To calculate Gy , we sum up the contributions [e a ,h a ], [e u ,h u ] from all carriers and ions 
in P^ 2 ) — AP( 2 ) and add the contributions [e u ,h„] from all particles in AP( 2 ). We do not take 
into account the contributions [e ,h ] from AP^ 2 ), because they are used to calculate the field 

Gg>. 

If k is sufficiently small, then in the corresponding points of observation 

G?> = G ( t ] = G T , (16) 
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since the properties and velocity distributions of particles are identical in — AP^\ and 
we may neglect the contributions [e ,h a ], [e u ,h u ] from particles in AP^ and the contributions 
[e u ,h u ] from particles in AP( 2 \ as the bounds of all these contributions do not depend on k. 
The equation (JTfij) being satisfied, we shall need no more information concerning the fields G^ 
to prove the formula (jSJ). However, we must calculate the fields G# . We have 

Gg> = 0, (17) 

because all contributions [e a ,h a ], [e u ,h u ] from all particles in P 1 ^ have been already used to 
calculate the field G^ . It remains to calculate G#\ 

3.3. The field G% ] 

Let points x' with coordinates x[, x' 2 , x' 3 in the Cartesian system X denote points belonging 
to p( 2 \ We define for a positive v that P^ is given by formulas: x[, x' 2 € [— v, u], x' 3 E 0], 
and that d 2 P^ lies on the plane x' 3 = 0. Let K be a unit vector normal to the plane £3 = 
and pointing in the positive direction of the a^-axis. 

Assume for the sake of simplicity that all electron carriers move only in the a^-direction, 
and that all of them have the same speed u outside the boundary layer AP^ (the results 
obtained by means of this assumption may be easily generalized for some more general cases). 
In the layer AP^ each electron stays during a time t K while it is accelerated from the velocity 
+«K to the velocity -mK. Integrating e (x, t) from (|12p with respect to time during the 
time t K , we obtain the averaged in a unit time contribution 5E^(x) to the field from one 
carrier. Let R = x, R — |R|, n = H/R (R is the distance from the center of the face 
&2P^ to the observation point). If R ^> v + k then with a good accuracy Ri = R and 
5Eg } (x) = -ZLxuR- 1 -!!*^ x K). 

If the number density of electron carriers in P^ is $o ; then the number of collisions 
with the layer AP^ per unit area and time is uQq/2. Therefore, we obtain E^(x) = 
— AS* Li& u 2 R~ 1 nx (n x K) where AS* is the area of d 2 P^ . Taking into account we 
can write 

Eg } (x) = -AS LET 1 nx(nxK) (18) 

where 

L = (2TTme c 2 )- 1 Q%W' J W' = mu 2 /2, (19) 

m is the mass of electron. 

Similarly, integrating h a (x, t) from (fT2|) . we get 

Hg } (x) = AS ceoLR- 1 nxK. (20) 

(2) 

Using the above method of derivation, one can also obtain the expressions for G^ in the general 
case when the electrons move in 3 space directions and have different velocities. In such a case, 
one gets the same expressions (fTHj) and (|2l]j). The only difference is that W in (fT§j) should be 
replaced by 

W = (1/2) kTrj(T) (21) 
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where W is the average kinetic energy of electrons per one degree of freedom, k is Boltzmann's 
constant, T denotes the temperature of electrons in kelvins, and n is a dimensionless coefficient 
equal to one for particles having the classical distribution of velocities. Some values of L and 
n for copper, silver, and a deuteron plasma are given in Appendix. 

Suppose that the values of $o an d u or W for electron carriers in cubes P^ are similar to 
those for carriers in copper at 300 K. Using the data from formulas (jA7|) in Appendix, we get 
then from ()18|) and (J2*U|) the numerical values of E#, mentioned in Section 2 for copper 
cube. Therefore, we obtain 

|Gg } |>0. (22) 



3.4. Conclusions 

From ©, JTSJ, JIBJ), and flI7J), we get 

G (1) = G T , G (2) = G r + Gg } . (23) 

Hence, taking into account (J22|) . we obtain that at least one of the following inequalities must 
be true: 

|G (1) |>0, |G (2) |>0, (24) 

which proves the formula 

As it has been already pointed out in Section 2, the result (|2"4"j) does not agree with the re- 
sults which can be obtained by means of the presently used methods of macroscopic averaging. 
These methods lead in our case to the conclusion that G« = G? = (i = 1, 2), because m 
both cases the macroscopic current and charge densities are equal to zero. One of the reasons of 
this contradiction is the fact that when we calculate the macroscopic field generated by carriers 
from PW, the contributions to this field from carriers leaving or entering Pw are discontinuous 
in time and do not satisfy the pertinent microscopic Maxwell equations. This invalidates, e.g., 
the indicated in Section 1 formulas given by de Groot (1969), and de Groot and Suttorp (1972). 



4. FIELD E BS GENERATED BY A SURFACE dP s 
IN AN INHOMOGENEOUS CONDUCTOR 

Let Ro , ip (Rq 3> ■0) be positive constants. Consider a cube P' defined by 

^e[-2Po,2Po], i = 1,2,3. (25) 

The cube contains an inhomogeneous alloy of silver and copper. On every surface x' 3 = 
constant the alloy is homogeneous. For x' 3 < —ip the medium in P' is a homogeneous copper 
while for x' 3 > ip the medium in P' is a homogeneous silver. Let 2-0 = 0.01cm, R = 1cm. 
Consider in P' a cylinder P5 defined by 

K) 2 + (4) 2 < Ro, 141 < 1-1^ (26) 
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We denote the lateral surface of Ps by 8lPs and the bases of Ps by db\Ps an d d^Ps ■ The 
base dbiPs lies on the surface x' 3 = —l.lip. Assume that the boundary of Ps is impenetrable 
for electrons. We shall calculate the boundary acceleration field E^s generated by electrons 
from P (P = P' — Ps) accelerated at the boundary of Ps- Let the observation point x lie in 
the center of Ps, i.e. x = (0,0,0). Let Ebsi, ^BS2, ^bsl be the contributions to E BS 
generated by electrons accelerated at dbiPs, d^Ps-, an d diPs, respectively. 

Integrating formula (JTHj) with respect to S over the base db\Ps in suitable polar coordinates, 
one obtains: 

E BS1 = 2nLKR m \R 01 - l.l^) 2 where P m = (R 2 + l.lV) (1/2) (27) 

(a more detailed derivation of formulas ([27]) to (|30|) is given in sections 5.2 and 5.3 of my 
previous pertinent paper, available electronically). In our case (Rq if)) we get with a good 
accuracy E BS i = 27rL Ctt _R K where L Cu is the value of L for electrons in copper. Similarly, 
E BS 2 = -2nL Ag R ~K. Hence 

Essb = E B 5i + E B S2 = i-2-RoK, L 2 = 2tt{L Cu — La 9 ) ■ (28) 

If k is sufficiently small, then it can be easily proved that the field E B sl can be estimated 
by integrating formula (|18|) with suitably changed K over diPs- In our case (2if) = 0.01 cm, 
R = 1 cm) one can obtain an estimate |Es5i|/|EB5b| < 3xl0~ 4 if one puts $ m the integrand 
twice as big as $ f° r copper. Therefore, with a good accuracy 

Ess = E BS b + E BS l = £ 2 -R()K. (29) 

From (J2HI) and (JM| we get L 2 = -3.9 x 10 15 Vm" 2 at 300 K and 

|E BS | > 3.5xl0 13 Vm -1 (30) 

for Rq = 1 cm . We see that the field E^s generated by particles from P' — Ps accelerated in 
the immediate vicinity of the whole boundary dPs may reach surprisingly great values even 
when the diameter of Ps is small. 

5. REMARKS 

This paper is a shorter and simpler version of my earlier paper described in Introduction. 
These of the more important results of that paper which have not been derived in the present 
paper are concisely formulated without proofs in the following remarks. Some methodological 
remarks have been also retained and are given in the next section. The earlier paper is available 
from the author on request electronically. 

Remark 1. It can be proved that formulas (|18j) . (J20j) with expressions (|19J1 . where W 
is replaced by W from (|21j). are also valid for a general medium where the thickness of the 
boundary layer k is comparable to, or greater than, the mean free path I of carriers in P (the 
definition of k in such a general case is given in the earlier paper). They are also valid for curved 
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boundaries when the moduli of radii of normal curvatures are sufficiently big in comparison 
with k and /, and in some cases when the latter condition is not satisfied. 

Remark 2. If the ions in P move and collide with the boundary c^P, then we must also 
take into account the boundary acceleration field produced by ions. This field is described by 
the same formulas (fTHjl. (f2Uj) with the data for electrons in expressions (fT^j) and (j2"Tj) replaced 
by the data for ions. Let the subscripts i and e label the quantities referring, respectively, to 
ions or nuclei and electrons. If, for instance, all ions have the same charge Qi and mass m i; we 
get from (JTHJ) and J2Q 

\Li/L e \ = (rji/rje) m e /mi 

as |Qj$oi| — IQe^oel- in this and in the other cases, the contribution to E B , H B from nuclei 
or ions may be usually neglected in comparison with the contribution from electrons, since the 
ratio m e /mj is small and r\ijr\ e < 1. 

Remark 3. Formulas (fTHj) . (J20J) with (|TT?|) and (}2"T]) result from expressions (fT2"|) . The 

latter are the first terms of the expansions with respect to the powers of u/c of the exact 
expressions. Therefore, (|18|) and (|2U|) are also the first terms of the expansions of the exact 
boundary acceleration field G^. If P contains a homogeneous medium, then the formulas (fT%|) . 
()20|) integrated over the closed surface dP vanish. In our cases the values of L are very large 
and the values of u/c are not very small (see Appendix). Therefore, one should consider the 
influence of further terms of the expansions of G^. Some partial results suggest that the second 
terms of these expansions integrated over dP may be different from zero for plasma in magnetic 
traps, even when the plasma were homogeneous. 

Remark 4. Let G ^ 0. Assume that G = G i + G 2 where G i does not depend on t 
and | G02 1 is sufficiently small. In general, the density J and the field Goi may influence the 
velocity distribution functions of particles in P. Hence the fields Gb and Gt may depend on 
J and p. However, if we assume that (@J) holds for J ^ 0, p ^ when &2P=dP or that |J|, \p\ 
are sufficiently small, then formulas Q and (fTUj) remain valid. 

6. METHODOLOGICAL REMARKS 

Remark Ml. The aim of this paper is a presentation of some new results deduced from 
the Maxwell equations with microscopic charge and current densities. These results concern 
averaged macroscopic electromagnetic field generated by large sets of moving charged parti- 
cles. We do not know the position and velocity of any particular particle, and may use only 
statistical information given by velocity distribution functions for various particle species in P. 
Nevertheless, we assume that the positions and velocities of particles are continuous functions 
of time. Such an assumption is necessary to use, e.g., the Lienard-Wiechert formulas. Without 
it any results concerning the electromagnetic field generated by sets of moving carriers are not 
a deduction from the Maxwell equations with microscopic sources even if they are an inference 
from some inductive extensions of these equations used in quantum theories. Therefore, we 
shall introduce a natural model based on this assumption and taking into account quantum 
considerations. 
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Let the region P contain N carriers with the same total electric charge Q and rest mass 
m, and let b denote a sufficiently small positive number. Let the macroscopic number density 
of these carriers be $(x, t) and their velocity distribution function be /(x, u, t). For sets of 
degenerate particles, such as electron plasma in metal conductors, the values of /(x, u, t) may 
result from quantum theories. We define in P a set of N privileged points with the following 
properties. The positions and velocities of these points are continuous functions of time, and 
their accelerations are piecewise continuous functions. Their trajectories do not cross d<iP ■ At a 
time to their macroscopic number density is equal to $(x, to)- If the distance between two priv- 
ileged points is less than 26, then they are strongly accelerated in random directions. Suitably 
choosing the distribution of such accelerations, we may ensure that the number density and 
velocity distribution of privileged points in P are given by $(x, t) and /(x, u, t), respectively. 

We regard each privileged point as a point belonging to the region occupied by a particle 
with the charge Q and rest mass m. When we write that particle p is at a point x at a time t, 
we mean that privileged point p is at x at t. All species of particles in P can be represented in 
this way. We may allow interactions between privileged points of different kinds of particles if 
we care about a more realistic description of collisions. 

Procedures equivalent to the setting of privileged points are universally used in statistical 
physics. For instance, they were explicitly or implicitly used by de Groot and Suttorp (1972), 
Ferziger and Kaper (1972), Landau and Lifshitz (1964), Schram (1991). 

Remark M2. Let d 3 x be a small volume in the system X in P, u n be the components 
of u in X, and d 3 u be a small volume in the velocity space u n . Assume that at a time t the 
expected number diV(x, u, t) of particles of a given species in the volume element d 3 x located 
at x, whose velocities lie in d 3 u about velocity u is given with a sufficient accuracy by the 
formula 

diV(x,u,t) = /(x,u) d 3 id 3 M 

where the velocity distribution function /(x, u) does not depend on t. 

Taking into account regions in the vicinity of the boundary c^-P, it may be sometimes 
convenient to consider subregions of P which are too small to contain a large number of particles 
(see, e.g., the next remark). We may then take into account a set of identical regions P n 
(n = 1,2,. . . J en ). In each region P m (m >1) the privileged points of all particles are located at 
a time t in the same points in which they were located in P m -\ at a time t- St. Considering all 
privileged points from all P m , we obtain for them the velocity distribution function f en equal 
to 7 en /(x, u). The reasoning may be then carried out for f en instead of /. The final result 
should be divided by I en . This approach makes dN well defined even when diV < 1. This is 
one of possible procedures making the concept of the velocity distribution function useful also 
in idealized or averaged microscopic regions. Such an averaging may be performed, e.g., in the 
microscopic vicinity of d^P on microscopic layers parallel to d<i P. 

Remark M3. We assume throughout this paper that the boundary d^P is impenetrable for 
carriers of electric charge. The form of the functions /(x, u) in a vicinity of d^P depends on the 
averaged distribution of forces acting on carriers in the neighborhood of c^-P. This distribution 
depends on the kind of medium in P and its theoretical idealization. We can also introduce 
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artificial boundaries and postulate extremal or not typical mechanisms of their impenetrability. 
In each case the forces acting in the vicinity of an impenetrable surface cause the averaged 
acceleration of particles to have a not vanishing resultant pointing into P. Consequently the 
number densities of carriers and the functions /(x, u) for various species of particles must, in 
general, depend on the distance from the point x to such a surface for x in a neighborhood 
of this surface. The only exception concerns the case when the barrier of forces making the 
boundary d^P impenetrable is sufficiently strong and steep to decelerate and accelerate the 
carriers within a distance which is very small in comparison with the average distance between 
particles in P. 

Remark M4. Consider a moving charged particle with the total electric charge Q. The 
electromagnetic field £ , H generated by this particle and observed at the point x and time t 
may be written as 

£(x, t) = e(x, t) + e'(x, t), W(x, t) = h(x, t) + h'(x, t), 

where the field e, h is equal to the field which would be generated in a vacuum by the charge Q 
located in the privileged point related to this particle. The field e', h' is a correction to the field 
e, h. This correction may be produced owing to the following reasons: the charge is not located 
in the privileged point, its density is finite, the charge density is strongly inhomogeneous and 
of opposite sign in some subregions of the particle, there are electric currents flowing and/or 
oscillating inside the particle, there are electric currents caused by the rotation of the whole 
particle. Let £ E (x, t), 7i E (x, t) be the sum of corrections e', h' from all particles in P, averaged 
in a short interval of time. The value of this sum depends on the kind of medium in P. In 
this paper we assume that if x is outside P and its immediate vicinity, then £ E (x, £) = 0, 
7i' E (x, t) = with a sufficient accuracy. This is equivalent to the assumption that outside P 
and its immediate vicinity we may neglect the total field of all electric and magnetic multipoles 
which may be introduced in order to describe the fields e', h' generated by all particles in P. 

Remark M5. The field Go generated by J and p in P is usually calculated by means of the 
retarded vector and scalar potentials A and tp. If, for instance, the medium in and outside P 
is homogeneous, and its permeability and dielectric constant are equal to 1, then the potential 
A is given by 

A(x, = — / d 3 x' |x - x'l- 1 J(x, t d ) 
AttJp 

where td = t — c _1 |x — x'|, and the potential tp(x.,t) is given by the analogical expression (see, 
e.g., Bochenek 1961). These potentials may be used to calculate the field Go if the Lorentz 
condition V • A + c~ 2 d(p/dt — is satisfied. This always takes place when the normal to dP 
component of J vanish identically on dP. It is a consequence of results given, e.g., by Stratton 
(1941) or Jones (1964). If this last assumption is not satisfied, then the field G is not defined 
while G remains well defined. 
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APPENDIX 

We shall calculate here the numerical values of L defined in (II 9 j) — with W replaced by 
W given in (j21J) — for electron carriers in copper and silver at a temperature of 300 K and for 
electrons in a deuteron plasma prepared for a controlled thermonuclear fusion reaction, with 
the number density of nuclei $od — 10 22 m -3 , the temperature of nuclei of 4x 10 8 K, and the 
temperature of electrons T = 10 7 K. The quantities referring to electrons in copper and silver 
will be denoted by subscripts Cu and Ag, respectively, while those referring to electrons in the 
deuteron plasma will have subscript D. 

The following constants are used: c = 3.00 x 10 8 ms _1 , ce = 2.65 x 10~ 3 Fs _1 , hp = 6.63 x 
10~ 34 Js (Planck's constant), k = 1, 38xl0~ 23 JKr 1 (Boltzmann's constant), m e = 9.11xl0 -31 kg 
(the rest mass of electron), = 6, 02x 10 23 mol -1 (Avogadro's number), Q = — 1.60x 10~ 19 C 
(the electric charge of electron) . 

We define 

W av = 3W= (3/2) kTi], (Al) 

(1/2) m e u 2 av = W av . (A2) 

We see that u av is an averaged speed of electrons and W av is their average kinetic energy in the 
non-relativistic range of velocities when these electrons move in all 3 space directions. It will 
be shown that for the electron plasma in copper and silver at 300 K we get c~ l u av < 5x 10 -3 . 
Therefore we may neglect the relativistic corrections to the mass and kinetic energy of electrons. 
As far as electrons in the plasma at 10 7 K are concerned, we receive c~ l u av < 0.08. In this case 
we also neglect the relativistic corrections, since this may cause at most a few per cent error of 
L D . 

Assume that the Fermi-Dirac statistics correctly describes the energy and velocity distribu- 
tions of electrons in electron plasma. If we consider electron plasma contained in a cube with 
a volume v in physical space, then the volume of the elementary cell in the momentum space 
is hp/ v. Each cell contains at most 2 electrons. Hence, if Co denotes the maximum energy of 
electrons at OK, and W 0av is the average energy corresponding to ( , we get 

Co = 4(8m e )- 1 (3$o/7r) 2 / 3 , W 0av = 0.6 Co- (A3) 

If kT/( <C 1, then we may neglect the increase of the average energy caused by temperature 
and put 

W av = W 0av . (A4) 

If kT/( ^> 1, then the Fermi-Dirac statistics becomes that of Maxwell and we may set 77 = 1. 
The exact analysis (see, e.g., Weizel 1955) gives for kT/( Q much smaller than 1: 

W av = W 0av [l + (5tt 2 /12)(A;T/Co) 2 ] + O[(kT/( ) 4 }. 

If fcT/Co is much greater than 1 then from formulas in Weizel's book one may infer (see also, 
e.g., Tamm 1957) that rj = 1 with a good accuracy if 

$0 hp(27im e kTy i/2 < 0.01 or T > 1.3 x 10~ 13 $ 2 /3 . 
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Therefore, for $od — 10 22 m~ 3 and T = 10 7 K, we get t)d = 1- From (|A1|) and ()A2|) we receive 
then c~ x u av B < 0.08. 

We calculate now Co for copper and silver. The mass of 1 mole of Cu is 63.5 g and of Ag is 
108 g. Assume that the specific mass at 300 K of Cu is 8.9 Mgm -3 , and of Ag is 10.5 Mgm -3 . 
Hence the number densities of atoms of Cu and Ag are, respectively, 8.4 x 10 28 m~ 3 and 5.9 x 
10 28 m -3 . The average number of electrons given to the electron plasma by one atom at 300 K 
is, respectively, 1.3 and 1.03 (see, e.g., Szczeniowski 1955, Clemmow 1973). Therefore, at 300 K 

$ 0Cu = Hxl0 28 m- 3 , $ A 9 = 6.1xl0 28 m- 3 , (A5) 
(ocu = 1.3 x 10- 18 J, W 0avCu = 0.78 x 10 -18 J, 

(A6) 

CoAg = 0.88 x 10- 18 J, W 0avAg = 0.53 x lO" 18 J. 
Formulas (jA6|) specify the maximum and average energies which the electrons would have at 

OK, were the numbers densities of electrons at OK equal to $oc« and &o Ag given by (|A5j) . 
For T=300K we have kT/( 0Cu = 3.2 x 10~ 3 , kT/( QAg = 4.7 x 10 -3 . Therefore we may use 
(jA4|) . From the approximate formula u av = {^WQ^m^ 1 ) 1 ^ 2 we get c _1 M ai ,c u < 4.6 xlO -3 , 
c~ 1 u av Ag < 3.9 x 10 -3 . It results from the Fermi-Dirac statistics that for nearly all electrons in 
copper and silver c _1 n < 7x 10~ 3 at T = 300 K. 

Taking into account flHJ), flHj), (|A"5j) . (|A"Hjl . we obtain from (0 and (jH} that at 300 K 

L Cu = -1.0xl0 15 Vm" 2 , ceqLcu = -2.7x 10 12 Am~ 2 , 

(A7) 

L Ag = -3.8xl0 14 Vm" 2 , ce L Ag = -1.0 x 10 12 Am" 2 . 

On grounds of ()A1|) . (jA3|) . and ()A4|) we notice that L Cu and L Ag should not practically depend 
on temperature in an interval within which the number densities flA5|) remain nearly constant. 
Replacing W av in ()A1|) by Wo av from ()A6|) . we receive for T = 300 K 

r] C u = 126, r]Ag = 85. (A8) 

Taking into account r\u = 1, we obtain for T = 10 7 K 

L D = -2.4xl0 10 Vm- 2 , ce L D = -6xl0 7 Am -2 . (A9) 
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